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Abstract. Very high energy gamma-ray data obtained by CAT and HEGRA from active galactic nucleus Mkn 501 
are used to constrain the cosmic Mid-Infrared background. While the entire infrared and submillimeter spectrum 
shape based on models has been fixed and the density scaled as a whole in previous studies, recent measures on 
the low and high energy infrared background are extensively used in this paper. In this original approach the 
infrared distribution is only varied in the unexplored 3.5-100 /an region. With conservative hypothesis on the 
intrinsic spectrum of Mkn 501, an upper limit of 4.7 nW.m _2 .sr _1 between 5 and 15 /im is derived, which is very 
close to the lower limit inferred from deep ISOCAM cosmological surveys at 15 /»m. This result is shown to be 
independent of the exact density of the A < 3.5 /im and A > 100 /im infrared distribution within the uncertainties 
of the measurements. Moreover, the study presented here rules out a com plete extragalactic origin for the 60 
micron tentative background detection as found by Finkbeiner et al. (200C). 
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1. Introduction 

Our knowledge of the early epochs of galaxies has recently 
increased thanks to the observational evidences provided 
by UV/Vis/Near-IR, far-IR and submillimeter (submm) 
surveys of high-redshift objects. In a consistent picture, 
galaxy formation and evolution can also be constrained 
by the background radiation which is produced by the 
line-of-sight accumulation of all extragalactic sources. 

In the last two years the Extragalactic Background 
(EB) at visible, IR and submm wavelengths has been 
finally constrained by both very deep source counts and 
upper limits on the diffuse isotropic emission at shorter 
wavelengths, and measurements in the submm range (see 
Gispert et al. 2000 for a review). 

The Cosmic Infrared Background (CIB) detected in the 
COBE data at wavelength greater than 100 /im (Fixsen et 
al. |1998|, Ha user et al. |1998| , Schlegel et al. |1998| , Lagache 
et al. 1999) contains a surprisingly large fraction of the 
background due to distant galaxies. The spectrum implies 
that the submm part of the CIB cannot be dominated 



by the emission of the galaxies which accounts for most 
of the CIB at 150 /im, and thus contains a unique 
information about high redshift IR galaxies. Considering 
the variety of long wavelength spectra observed for 
these galaxies, Gispert et al. (200C) have shown that a 
co-moving production rate of far-IR radiation with strong 
evolution at low redshifts but little evolution between 
redshifts 1 and 4 is the only solution allowed by the CIB. 
The search for the UV, optical, near-IR and mid-IR 
EB (from 0.2 to 2.2 /mi) obtained by summing up the 
contributions of galaxies using number counts currently 
gives only lower limits; but the flattening of the faint 
counts suggests that we are now close to convergence. 
Nevertheless, number counts give a good determination 
of the background only if the full flux of the sour ces is 
properly taken into account (Bernstein et al. 1999] ) and 
if diffuse components are negligible. In fac t, stu dies of 
the diffuse COB E emission (Gorji an at al. [2000 ; Dwek 
& Arendt 1998 ; Cambresy et al. 2000 ) suggest that a 
significant fraction of the EB is missed from number 
counts integrals in the near-IR. 
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Due to the bright zo diacal foreground emission (see 
for example Kellsall et al. 1998 ), EB determinations from 
5 to about 80 /tm are very difficult. In this wavelength 
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range, EB constraints come mostly from deep IR cosmo- 
logical surveys as the ones performed with ISOCAM (e.g. 
Elbaz et al. 1998). A more controversial tentative CIB 



detection at 60 /im has also been published using DIRBE 
data (Finkbeiner et al. |2000). 



It has been suggested that observations of the TeV 
spectrum of extragalactic sources can be a powerful tool 
to constrain the CIB spectrum (Nikishov 1962; Gould 
& Schreder 1967; Stecker et al. 1992), especially around 
10 /im, where the EB constraints are today very weak. 
As a matter of fact, TeV gamma rays propagating in 
the intergalactic medium undergo absorption through 
electron pair production on CIB photons. Using sources 
Mkn 421 and Mkn 501, meaningful u pper limits have 
been established (see Malka n fe S tecker 1998 , Funk et al. 
1998, Macminn fc Primack 1996, Stanev & Franceschini 



1998, Guy et al. 2000 for example) 



In this paper, gamma-ray data are used, together 
with the best presently available direct CIB measurements 
(rather than CIB models), to improve the constraints in 
the most difficult spectral region, near the maximum of the 
zodiacal emission. The paper is organised as follows. We 
discuss in Sect. 2 the CIB measurements that are going to 
be used. Sect. 3 presents the TeV gamma-ray data. The 
determination of the upper limit is explained and done 
in Sect. 4 and 5. In Sect. 6, we compare our upper limit 
with lower limits derived from deep source counts obtained 
with the ISOCAM instrument aboard ISO (Elbaz et al. 



1998, Metcalfe et al. 2000) and discuss the CIB distribu- 



tion shape from 10 to 100 /im. Cosmological implications 
are then discussed. Finally, in section 7, the consequences 
on the physical parameters of Mkn501 regarding the in- 
ternal absorption are briefly considered. 

2. CIB density distribution 

This study aims at deriving upper limits on the minimum 
of the CIB distribution lying in the Mid-IR (5-80 /im) re- 
gion. In the Far-IR region, the distribution is constrained 
to account for the 100 /im determination from Lagache 
et al ( |l999| ), which is in ve ry go od agreement with the 
value of Finkbeiner et al. (2000); data at longer wave- 
lengths are not relevant for this study as they are be- 
yond the threshold of the hardest TeV photons. In the 
Near-IR region, three different hypothesis, parametrised 
in this work by polynomial fits between 0.2 and 3.5 /im, 
are considered: the first one (dashed-dotted line in Fig. |l|, 
referred as HDF) is directly based on Hubble Deep Field 
galaxy counts (Pozzetti et al 1998), and can therefore 
be considered as a lower limit; the second one (full line 
in Fig. [l], referred as HDF + spectro) is derived from 
HDF results co mbined with ground-based spectrometry 
(Bernstein et al. 1999|) and the third one (dotted line in 
Fi g. |l], r eferred as Cambresy) is based on (Cambresy et 
al. 2000] ) points at 1.25 and 2.2 /i m linked to the UV point 
at 0.2 (im (Bowyer et al. ( 1991 ) combined with Armand 



et al. (1994)), as it is clearly incompatible with HDF and 
HDF + spectro hypothesis^. All these parametrisations 
are required to agree with the 3.5 /jm measurement from 
Gorjian et al. fl2000| ); HDF and HDF + spectro hypoth- 
esis also include the 1.25 and 2.2 /im points from Wright 
( p000| ). 

The CIB coming from the line of sight accumulation 
of all extragalactic sources is expected to have a smooth 
spectrum. The CIB might present a minimum between the 
stellar dominated part and the interstellar dust dominated 
one between 3 and 15 /im. To make optimum use of the 
upper limits that TeV gamma rays can provide we thus use 
an empirical model for the CIB with a small number of free 
parameters (X m in and f m in) which interpolates between 
3.5 and 100 /im two segments of parabola (in log-log plot) 
linking the 3.5 or 100 /im points to the minimum located 
at (Xmin, fmin) with a null first derivative at this point. 
The CIB empirical density distribution is shown in Fig. ^ 



for A r . 



10 /im and /„ 



2, 4 or 6 nW m 2 1 



3. TeV gamma-ray astronomy data 

The TeV source Mkn 501 is the second closest X-ray 
selected BL Lac object after Mkn 421 with a redshift 
z ~ 0.034. During the 1997 outburst which lasted sev- 
eral months, Mkn 501 was ob served intensively in X- 
rays (B eppoS AX: Pian et al. |l998[ RXTE: Lamer & 
Wagn er 1998 ) and TeV 7-ray s (Whipple: Catanese et al. 



1997 , Samuelson e t al. |1998| ; HEGRA: Aharonian et al. 
1997 , |!999a| , |l999fc| ; Telescope Ar ray: Hayashida et al. 
1998; CAT: Djannati-Ata'i ct al. 



1999[ ). The exceptional 



April 16, 1997 flare was observed by BeppoSAX and low- 
energy threshold (~ 300 GeV) Whipple and CAT atmo- 
spheric Cherenkov telescopes, allowing the derivation of 
the energy spectrum with good accuracy in a broad dy- 
namical range. 

These unique data initiated interesting efforts to set 
mean ingful upper limits on th e CIB flux (se e Bill er et al. 
1998|, Stanev & Francesch ini |l99S| , Barrau |l998| , Stecker 



& De Jagcr |1998| , St ecker [1999| , Aharonian et al. |l999b 
Coppi fc A haronian 1999 , Konopelko et al. 1999 , Guy 
et al. |200C| ). 

In this paper, CAT and HEGRA data from 400 GeV 
to 17 TeV are used. Due to bad weather, the HEGRA 
telescope system could not observe this flare. But the 
HEGRA observations of Mkn 501 revealed that, despite 
flux variations in sub-day time-scales, the shape of the 
energy spectrum above 1 TeV remained essentially sta- 
ble throughout the entire state of high activity in 1997. 
The H EGRA 'time-averaged' spectrum (Aharonian et al. 



1999b[ ) can therefore be used after normalisation to the 
CAT April 16, 1997 spectrum at E =1 TeV with the re- 



scaling factor of «2.2 (Guy et al. |2000|) . At these energies 
the agreement between the CAT and HEGRA spectra is 
quite impressive. Below 1 TeV, CAT points are used solely 



1 However, we will see (Sect. 5.2) that the impact of the 
density below 1 /im is «null 
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Fig. 1. Assumed shapes of the CIB density distribution 
between 0.2 and 100 /im. Each profile between 3.5 and 
100 /im is characterised by the position of the mini- 
mum (Xmin, fmin)- I n this example A m i„ = 10 /im. See the 
text for details. Arrows show the 90 % interacting inter- 
val of target radiations for a 4 or 17 TeV photon (for 
fmin= 4 nW m~ 2 sr _1 ). The full line links together the 



(HDF+spectro) points (0.3 to 0.8 /im), the (Wright 2000 ) 
points (1.25 to 2.2 /im) and the (Gorjian et al. 2000 ) 
points (2.2 to 3.5 /im). The dashed-dotted line lin ks to- 
gether smoothly t he HD F points to the (Wright 2000 ) 
and (Gorjian et al. 2000[ ) points. The dotted line links to- 
gether the UV measurement at 0.2 /im, the (Cambresy et 
al. |2000|) p oints (1.25 to 2.2 /im) and the 3.5 /im (Gorjian 
et al. 2000) points. 



since in this energy region both the statistical and system- 
atic errors of the data obtained close to the energy thresh- 
old of the HEGRA telescope system are significantly larger 
than CAT data uncertainties. 



4. Unfolding TeV spectra 

The influence of low energy photons in the Universe on 
the propagation of Very High Energy (VHE) gamma-rays 
was pointed out by Nikishov ( 1962 ). An original way of 
using ground-based TeV observations of dist ant so urces 
to probe the CIB was given by Stecker et al. (1992). The 
fundamental idea is to look for absorption in the intrinsic 
spectrum as a result of electron-positron pair production 
by photon collisions jTeV + Icib — * e + + e~. 

In such an interaction between a gamma-ray of 
energy (1 + z)E and an infrared photon of energy 
(1 + z)e, where z is the redshift, E and e the observed 
energies at 2 = 0, the pair production threshold is 



between photons and m the rest mass of the electron 
The cross-section can be written as (Heitler 1960) 



= k(l - 



P 2 ) (2/3 (f3 2 



-2) + (3-f3*)ln i±f cm 



1/2 



with P= (1-2 {mc 2 ) /{Et (1 - cos8)) (1 + z) J and 

k = 1.25 x 10~ 25 . If the infrared photons have a density 
number n (e) de cm -3 , the corresponding optical depth 
for attenuation is 



t{E) 



B£ So dz (! 
J_. d (cos9) 



z) 



1/2 



1 — COS0 



I™ de n{e)a{E, e, 



where e t = 2{mc 2 ) 2 /{E{1 - cos6){\ + z) 2 ), z s is the red- 
shift of the source, c the light speed, and H the Hubble 
parameter (assumed here equal to 65 km s _1 Mpc _1 ). For 
this equation, the density parameter fin. was chosen close 
to unity. The detected flux is then attenuated by a fac- 
tor e~ T ( E \ The CIB energy distribution is assumed to be 
independent of z as the 7-ray source redshift is very low 
(0.034). 

The maximum cross-section is reached for an infrared 



photon wavelength of Xci b ~ A c 



2mc- 



where X r 



h/ {mc) 



Ee{l + zf (1- cosB) > 2( 1 



where 9 is the angle 



is the Compton wavelength of the electron. As shown in 
Fig. [j], 7-photons with energy between a few TeV and 
20 TeV "see" CIB photons with wavelengths between 3.5 
and 100 /im. 

5. Upper limit on the CIB 

5.1. Method 

To turn out the absorption corrected spectra into an 
upper limit on the CIB density, the intrinsic spectral 
energy density of Mkn 501 is assumed to be convex in 
the multi-TeV region (Guy et al. 2000 ). This conservative 
hypothesis is based on the fact that no natural physical 
process can re-inject energy above the Inverse-Compton 
bump maximum. Both Klein-Nishina effect and auto- 
absorption within the source would cause the flux to 
decrease more and more rapidly as a function of energy. 
In the framework of more complex hadronic models this 
convex shape is also expected, either because of proton- 
initiated-cascade (often modelized as a broken power-law 
around 3 TeV), or due to the inclusion of /i-synchrotron 
radiation (Mannheim, private communication). Even in 
the extreme case where proton synchrotron radiation is 
at the origin of the TeV bump (Aharonian |2000 ), the 
shape of the emission remains convex. 

Looking at the absorption corrected spectra shown 
in Fig. |^, it cannot be excluded that the maximum of 
the so-called Inverse-Compton peak is not yet reached 
at 17 TeV when a large amount of CIB radiation 
{fmin > 6 nW m~ 2 sr _1 ) is assumed. Physical parameters 
(essentially magnetic field and Doppler factor) required 
to produce such a spectral energy distribution maximum 
above 17 TeV are substantially disfavoured but, even 
in this case, the vF{v) shape should remain convex. 
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Fig. 2. Shapes of the reconstructed TeV spectra. The open 
circles show the observed CAT+HEGRA TeV spectrum 
while the filled symbols show the absorption corrected 
spectra (HDF + spectro hypothesis). The results of the 
parabola fit are superimposed. 



This latter point is ensured by B eppoSAX 1-100 keV 
measurements (Pian et al. 1998) showing a clearly 
convex spectrum before the synchrotron bump maximum 
(around 100 keV) which is supposed to be mimicked by 
the TeV-spectrum before the so-called Inverse- Compton 
bump maximum. This hypothesis is reinforced by the 
fact that the sub- TeV slope, which is independent of the 
CIB density beyond 3.5 /tm, effectively reflects the X-ray 
slope in the keV range. This behaviour (together with the 
correlated variability reported e.g. by Aharonian et al. 
( |1999a| ) and Djannati-Atai et al. ( |1999| )) indicates that 
the same population of particles is at the origin of both 
X-ray and 7-ray emissions, whatever this population is. 
In particular, the self-synchro-Compton model fi ts sat - 
isfactorily the absorption corrected data (Guy et al. 2000| ) . 



In order to quantify the concavity of the absorption 
corrected spectrum, a parabola fit is performed from 
3.3 to 17 TeV in the plane (log(v),log{vFv)). This 
function, which is simple, is chosen for its constant second 
derivative a = d 2 {log(vFv)) / 'd 2 (logiy)) , avoids the choice 
of a particular test-energy and fits satisfactorily the data. 
The previous physical constraint on the concavity of the 
TeV spectrum simply reads as a < 0. The parameter a 
is computed for f min in the range 0.1-7 nW raT 2 s _1 . 
Fig. U shows the result of the parabolic fit in several cases, 
superimposed to the absorption corrected experimental 
points. The fit under-estimates the energy density of the 
hardest photons, ensuring that the local second derivative 



is under-estimated and thus making the test conservative. 

Furthermore, when a > 0, fitted parabola with 
-b/2a > 17 TeV, where b = d(log{vFv))/d(log{v)), are 
accepted. This ensures that, whatever the concavity, no 
fitted function decreasing with v over the energy range 
considered here will be excluded. Asking for a decreasing 
vF(v) distribution is the most natural constraint in the 
previously mentioned (Inverse-Compton bump maximum 
below or around 1.4 TeV) most favoured case. This double 
test on both the first and the second derivative makes the 
resulting upper limit very conservative. 

5.2. Reliability of the upper limit 

Fig. |^ shows the impact of the main uncertainties on the 
measured part of the CIB density distribution on the ab- 
sorption corrected TeV spectrum. The observed spectrum 
of Mkn 501, combining CAT and HEGRA data, is plot- 
ted with open circles in panel a). Spectra corrected from 
CIB absorption are superimposed with filled symbols in 
different cases: a) HDF, HDF + spectro and Cambresy 
hypothesis b) 3.5 (im ± 1 a (NIR); c) position of the min- 
imum \ n in = 5 or 15 /xm; d) 100 /im value ± 1 a (FIR). 
By default, A m ;„ = 10 /im. 

From panel a), it can be seen that the impact of 
the NIR density below 1 /tm is null beyond a few TeV, 
as the threshold for pair production for a 1 TeV pho- 
ton is approximately 1 /im. Furthermore, even below 1 
TeV, HDF and HDF + spectro hypothesis lead to the 
same absorption corrected spectrum. Therefore, only the 
HDF + spectro and Cambresy hypothesis will be con- 
sidered hereafter. The value of \ m in (panel c) and the 
FIR density at 100 /im - in a minor way - (panel d) can 
slightly change the very high energy tail of the source spec- 
trum (above 10 TeV) while uncertainties on the 3.5 /tm 
measurement have no sizeable impact (panel b). So, the 
uncertainty on the 100 /tm measurement and a range of 
values for the parameter X m i n (5-15 /tm) must be taken 
into account when deriving upper limits on the MIR den- 
sity. 

The evolution of the parameter a in function of f m in 
is presented in Fig. |] in the case X m in — 10 /tni and a 
HDF + spectro CIB density distribution while the values 
of fmin corresponding to a null parameter a (within 5a) 
for several values of X m in are gathered in Tab. [l] for all 
hypotheses regarding the NIR or the FIR density. 



5.3. Conclusion 

The results of this work are given in Tab. § The 3.5- 
100 /tm photons being essentially absorbed by 2-15 TeV 
gamma-rays, the CAT and HEGRA spectra of Mkn 501 
are particularly well suited to test the MIR density. The 
high flux allows a good spectral resolution and the simul- 
taneous X-ray data provides crucial information on the 
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Fig. 3. Impact of the shape of the CIB spectrum on the 
absorption corrected TeV spectrum for i m in= 5 nW m~ 2 
sr -1 . The open circles (panel a) show the observed TeV 
spectrum while the other symbols show the absorption 
corrected spectra with: a) HDF, HDF + spectro and 
Cambresy hypothesis (resp. stars, filled circles and open 
crosses); b) 3.5 /mi value - or + 1 a (resp. stars and filled 
circles); c) position of the minimum X m i n = 5 or 15 /im 
(resp. stars and filled circles); d) 100 value - or + 1 a 
(resp. stars and filled circles). By default, X m i n — 10 /im. 
Error bars are hidden by the experimental points. 
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Fig. 4. a parameter of the parabola fit versus f m in 
(log{vF(v)) = a x log 2 {v) + b x log(v) + c) with 5 a errors 



(for \ m in — 10 /im and a HDF 
distribution.). 



spectro CIB density 



in the range 5-15 /im, the Mid-IR density upper limit re- 
mains lower than 4.7 nW m~ 2 sr -1 . 

Our results are in agreement with previous works: for 
example, Stanev & Franceschini (1998|) obtained an up- 



per limit of 4 to 5 nW m sr at 15 /tm with two 
CIB hypotheses and Stecker & 
an upper limit of 4 nW m~ 2 sr 



De Jager (1997) derived 
1 at 20 /im. It shows the 
" robustness" of the method based on TeV gamma-rays ab- 
sorption by IR photons, essentially due to the peak in the 
pair production cross section. The NIR and FIR fluxes are 
now confirmed by direct measurements and these obser- 
vational results enhance the reliability of the constraints 
on the MIR intensity distribution. 

An approach similar to ours has also been proposed by 
Dwek et al. (2000) who found a CIB density lower than 



5 nW m 2 sr 1 in the range 6-30 /tm and lower than 
10 nW m -2 sr -1 at 60 /im. 



synchrotron bump shape. The MIR upper limit derived in 
this paper is robust as it is alomost independant of the 
NIR hypothesis (within 1 a error), and only marginally 
depends on the FIR flux hypothesis (within 1 a error). 
Somewhat more surprising is that the exact position of 
the minimum in the 5-15 /tm range significantly affects 
the upper limit (as seen in Fig. ||). The quite high sen- 
sitivity to the value of X min is due to the fact that the 
shape of the 7-ray spectrum strongly depends on the CIB 
density distribution between 5 and 50 /tm (and not only 
on its mean intensity). In other words, more absorption is 
allowed if the mid-IR spectral distribution, and thus the 
energy dependence of the extinction, gives an acceptable 
shape of the TeV spectrum of the source. 

This limit is conservative, both due to the method 
which tends to under-estimate the high energy flux density 
and to the weak hypothesis regarding the spectrum con- 
cavity. Indeed, internal effects, such as Klein-Nishina cut- 
off and self-absorption, lead to very sharp multi-TeV in- 
trinsic spectra. The upper limit is more constraining than 
the indirect upper limit derived from the same data on 
Mkn 501, which relics on a global CIB scaling (Guy et al. 
200C ). Moreover, whatever the position of the minimum 



Table 1. 5 a upper limit on the Mid-IR density as a 
function of the position of the minimum X m i n and of the 
CIB density distribution. 



IR density 


5 a Mid-IR density upper limit in nW m M sr 1 
^min~5 /tm 10 /tm 15 /tm 


HDF + spectro 


1.3 


2.8 


4.2 


Cambresy 


1.2 


2.6 


3.9 


100 /tm -1 a 


1.8 


3.4 


4.7 


100 ^m +1 a 


1.1 


2.5 


3.9 



6. Discussion on the CIB measurements 

6.1. The case of the 60 micron tentative EB detection 

From analysis of t he DI RBE weekly averaged sky maps, 
Finkbeiner et al. ( 200C ) have detected substantial flux 
in the 60 /tm DIRBE channel in excess of expected 
zodiacal and Galactic emission. The observed signal is 
consistent with an isotropic background at the level 
vl v = 28.1 ± 1.8 ± 7(syst) nW m~ 2 sr" 1 . While this new 
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excess is not necessarily the CIB, they have ruled out all 
known sources of emission in the solar system and Galaxy. 
They therefore tentatively interpret this signal as the 
CIB. However, they point out that the IR excess exceeds 
limits on the EB derived from the inferred opacity of the 
inter-galatic medium to observed TeV photons. 

In the assumption that the residual observed at 
60 /im is the CIB, we can compute, using the same 
method as previously, a Mid-IR upper limit on the 
CIB. With A m i„ = 10 /xm, we obtain an upper limit of 
1.0 nW m~ 2 sr _1 which is clearly incompatible with the 
lower limit of 3.3 nW m~ 2 sr" 1 at 15 /im, derived from 



deep ISOCAM surveys (Biviano et al. 200C). Even when 



the 60 /im tentative measurement is decreased by 1 c, 
the absorption corrected TeV spectrum in unacceptable 
(upper limit of 1.3 nW m~ 2 sr _1 ). With A,„m=15 /im, 
the upper limit becomes marginally compatible with the 
ISOCAM result (upper limit of 2.2 nWm" 2 sr" 1 ). Thus, 
as pointed out by Finkbeiner et al. ( [2000 ), we conclude 
that the 60 /im excess is certainly not only CIB. It con- 
tains residual foreground emission, more likely zodiacal 
residual emission (see Sect. 6.3). 



6.2. The EB around 15 \im 

Fig. ^| shows together our results with the 6.5 and 15 /im 
lower limits obtained by integrating the energy from 
sources observed in ISOCAM deep fields. The flattening 



of the faint counts at 15 /im (Metcalfe et al. 2000) sug- 
gests that we are now close to convergence and thus that 
number counts are not very far from the true value of 
the EB (if no diffuse emission contributes to it) . Number 
counts give 3.3 ±1.3 nW m -2 sr~ - 1 (Biviano et al. p000| ). 
Using the model of galaxy evolution from Dole et al. (in 
prep), we find an EB of about 4.4 nW m -2 sr" 1 , which 
is very close to the upper limit derived using the gamma- 
rays for a high value of X m in- This tends to show that 
the absorption in the gamma-ray spectrum of Mkn 501 
is mainly, or entirely due to absorption through electron 
pair production on CIB photons (see Sect. 7). Thus, the 
upper limit obtained with A m i„ = 15 /im is probably a best 
guest of the CIB level. Using a different argument based 
on the expected intrinsic shape of the TeV spectrum and 
the "high" CIB model from Stecker & De Jager ( |199SD , 
Konopelko et al. ( 1999] ) also concluded that the observed 
curvature of the spectrum of Mkn 501 could be naturally 
explained by the extragalactic absorption. 



6.3. CIB shape between 10 and 100 fim 

In the previous sections, the CIB density distribution was 
assumed to be described by a parabola betw een the min- 
imum and the 100 /im Lagache et al. ( |l999| ) point. This 
function may underestimate of the EB light and it is im- 
portant to study the consequences of a change in the CIB 
shape. As illustrated in Fig. [g, two more realistic (and less 



10 




10 

A (/im) 

Fig. 5. CIB density in function of the wavelength: data 
and upper limits from VHE observation (HDF + spectro 
case). Triangles linked by the wide full grey line indicate 
the data used in this paper to constrain the NIR and FIR 
parts: below 3.5 /im the combination of HDF measure- 
ments and ground-based spectrometry (Bernstein et al. 
1999| ), (Wright |2000| ) and (Gorjian et al. ^000[ ) points 
and, above 100 /im, DIRBE results (Lagache et al. 
|1999D - The nu mber c ounts results from ISOCAM 15 /im 
(Biviano et al. |2000| ) and 6.5 /im deep cosmological ob- 
servations (Desert, private communication) are shown by 
stars while the squares come from DIRBE/COBE obser- 
vations (Hauser et al. 1998). The two other hypothesis 



below 3.5 /im (dotted lines) lead roughly to the same up- 
per limits. 



conservative) distributions have been tested with a mini- 
mum at (10 /im, 2.8 nW m -2 sr _1 ), i.e. in agreement with 
ISOCAM measurements. In each case, even taking into ac- 
count only the 140 /im Lagache et al. point, the shape of 
the absorption corrected TeV spectrum is rejected by the 
physical criteria given in section 5.2. This means that, at 



least, the 100 /tm Lagache et al. (1999) measurement may 
be overestimated. 

The DIRB E zod iacal emission model was obtained by 
Kelsall et al. ( |l998|) relying on its time variability^ This 
model is critical to get the CIB in the near infrared. Its 
accuracy can be estimated using the residuals observed at 
wavelengths were the zodiacal emission is maximum (12 
and 25 /im). The residual emission, obtained by Hauser 
et al. ( |1998| ), has in fact a spectrum very similar to the 
zodiacal one. The residuals are about 470 nW m -2 sr _1 at 
12 and 25 /im. The work presented just before reveals that 



2 Oth er mo re empirical models h ave also been built by 
Wright fl2000|) and Finkbeiner et al. (|2000|). 
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a) C II B density distributions 
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b) TeV unfolded spectra 
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Fig. 6. Consequences of a change in the CIB shape be- 
yond lO/im: hypothesis on the CIB density distribution 
(panel a) and corresponding absorption corrected high en- 
ergy Mkn 501 spectra (panel b). Ordered by increasing 
densities, we consider the normal {i.e. parabolic) distribu- 
tion (dots) and two cubic distribution linking the 10 /im 
minimum to the 140 /im Lagache et al point (squares) and 
the 100 /im Lagache et al point (up-triangles). 



this best model may underestimate the zodiacal emission 
as the uncertainties on other contributions (instrumental 
and interstellar) are significantly smaller. We can thus 
make a conservative estimate of the amount of zodiacal 
emission which is not removed by this model at 12 and 
25 ^m to be about 400 nW m~ 2 sr -1 . The amount not 
removed at 60 and 100 /im is thus 40 nW m~ 2 sr _1 and 
8.4 nW mT 2 sr -1 (using the Kelsall et al. 1998 smooth 



high latitude zodiacal cloud colour ratios). This reduces 
the CIB at 100 /im from 23.4 to 15 nW m" 2 sr" 1 . 
At 60 /im the extra-zodiacal emission to be removed 
(40 nW m -2 sr -1 ) is comparable to the residuals and 
thus no meaningful value can be obtained on the CIB 
at this wavelength. Using the same arguments, the 
extra-zodiacal emission to be removed at 1.25 /im in the 
Cambresy et al. ( |2000| ) work is 38.0 nW m~ 2 sr" 1 from 
the published result of 56.9 nW m~ 2 sr -1 . 

The impact of a larger zodiacal emission on the upper 
limits in the MIR range is shown in Fig. |^. With a flatter 
CIB distribution, the upper limit at 10 /im increases from 
2.8 to 3.7 nW m~ 2 sr -1 . It remains very close to the direct 
measurement from ISOCAM at 15 /im. 

In addition to a high interest for cosmology and galaxy 
evolution, a better knowledge of the entire CIB shape 




HDF+spectro 

+ Wright +Gorjion 



Fig. 7. CIB density vs wavelength: data and upper lim- 
its from VHE observation (A TO j n = 10 /im case). The 
HDF + spectro spectrum is linked to different points at 
100 /tm: Lagache et al. ( |1999[ ) flux (23.4 nW m~ 2 sr" 1 , full 
line), the same Lagache et al. flux with a higher zodiacal 
background subtracted (15 nW m~ 2 sr -1 , dashed line) 
and this point -1 a (8.7 nW m~ 2 sr -1 , dotted line). The 
number counts result from ISOCAM 15 /im (Biviano et al. 
[2000 ) is shown by the star. 



would be of great interest to unfold TeV spectra of ex- 
tragalactic sources. 

7. Discussion on Mkn501 physical parameters 

In the previous section, so as to derive a conservative up- 
per limit on the CIB density, no internal absorption within 
the source has been taken into account. Nevertheless, the 
ISOCAM 15 /im point (at 3.3 nW m -2 sr -1 ) gives a lower 
limit on the CIB density and, assuming this value to be a 
measurement, the resulting maximum possible absorption 
within the source can be estimated. The previous physical 
constraints on the absorption corrected spectrum can 
be applied again, assuming the minimum possible CIB 
density (to remain conservative), and varying the optical 
depth tagn due to 7 — 7 pair production within the jet 
ofMkn501. 

To compute tagn in the TeV range, the low-energy 
spectrum (i.e. the energy distribution of target pho- 
tons) has to be known. We use the approximation de- 
rived by (Bednarek & Protheroe, 1999] ) using a fit to the 
BeppoSAX observations made during the April 15/16 flar- 
ing activity together with an indication from OSSE ob- 
servation made during the high state in 1997 April that 
the spectrum continued to approximately 500 keV with 
the same energy flux per log energy interval. Performing 
a numerical integration of the optical depth formula (cf 
Sect. 4), the authors show that, in the framework of ho- 
mogeneous SSC model, : 

tagn « 3 x 10 s D- 4s E 0A t^ r 

where D is the Doppler factor of the blob where rela- 
tivistic electrons are assumed to be confined, t var is the 
variability time-scale of the source in seconds and E is the 
7-rays energy in TeV. The efficiency usually associated 
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Fig. 8. Optical depth of the TeV flux in function of 
the energy E. The stars symbolise the absorption due 
to the EB photons while the line shows the absorp- 
tion due to IR photons in the inner jet (A m i„=15 /zm, 
f mm =3.3 nW m- 2 sr" 1 ). 



with this mechanism has been supposed, as usual, equal 
to unity. 

Both absorptions are compared in Fig ||. They look 
rather similar as the dominant behaviour is in E~ 2 in en- 
ergy density. But the internal absorption is continuously 
increasing with energy while the CIB absorption exhibits a 
"plateau" around 1-4 TeV. Adding this new optical depth 
to the one associated with 7 — CIB interactions, we de- 



rive a 5cr upper limit for 3 x 10 D 



8 D-4.8+-1 



of 0.6 (com- 



puted with A m i„=15 M m an d /mm=3.3 nW mT 2 sr _1 , 
i.e. for the ISOCAM point). With a reasonable variability 
time-scale around 1 hour, the resulting minimum Doppler 
factor is D > 11.7. Despite the large CIB absorption and 
the resulting stringent upper limit on the internal opti- 
cal depth, these Doppler factors are in agreement with 
classical AGN descriptions and the simple homogeneous 
SSC model works satisfactorily. Nevertheless, it should be 
pointed out that this lower limit on the Doppler factor in- 
creases as the wavelength of the minimum of CIB density 
decreases and becomes infinite below w 10 /im when there 
is no space left for internal absorption, unless an inhomo- 
geneous jet is assumed. 

With the next generation of instruments (VERITAS 
and MAGIC in the northern hemisphere, HESS and 
CANGAROO in the southern one) which will start oper- 
ating in 2002/2003, good quality spectra of several AGNs 
at different redshifts, will be available. It will lead to, si- 
multaneously, a better understanding of the TeV emission 
mechanism and a better knowledge of the CIB density 
distribution. 
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8. Note added in proof 



In this paper an early version of the Cambresy et al. 2000 



paper was used: the 1.25 /im flux was overestimated by 
5 %: it does not affect our results. 



